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Abstract. The sidereal anisotropy of galactic cos- 
mic ray (GCR) intensity observed with the Tibet 
Air Shower (AS) experiment still awaits theoretical 
interpretation. The observed global feature of the 
anisotropy is well reproduced by a superposition 
of the bi-directional and uni-directional flows (BDF 
and UDF, respectively) of GCRs. If the orientation 
of the deduced BDF represents the orientation of 
the local interstellar magnetic field (LISMF), as 
indicated by best-fitting a model to the data, the 
UDF deviating from the BDF orientation implies a 
significant contribution from the streaming perpen- 
dicular to the LISMF. This perpendicular streaming 
is probably due to the drift anisotropy, because 
the contribution from the perpendicular diffusion is 
expected to be much smaller than the drift effect. 
The large amplitude deduced for the UDF indicates 
a large spatial gradient of the GCR density. We 
suggest that such a density gradient can be expected 
at the heliosphere sitting close to the boundary of 
the Local Interstellar Cloud (LIC), if the LIC is 
expanding. The spatial distribution of GCR density 
in the LIC reaches a stationary state because of the 
balance between the inward cross-field diffusion and 
the adiabatic cooling due to the expansion. We derive 
the steady-state distribution of GCR density in the 
LIC based on radial transport of GCRs in a spherical 
LIC expanding at a constant rate. By comparing 
the expected gradient with the observation by Tibet 
experiment, we estimate the perpendicular diffusion 
coefficient of multi-TeV GCRs in the local interstellar 
space. 
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I. Introduction 

The galactic cosmic rays (GCR) are high-energy 
nuclei (mostly protons) produced in our Galaxy. The 
directional anisotropy of GCR intensity gives us im- 
portant information on the magnetic structure of the 
region in space through which GCRs traveled to the 
Earth. In this paper, we analyze the sidereal anisotropy 
of ~ 5 TeV GCR intensity observed by the Tibet Air 
Shower (AS) experiment, which is currently the world's 
highest precision measurement of GCR intensity in this 



energy region, utilizing both the high count rate and 
good angular resolution of the incident direction. The 
GCR anisotropy in this energy region is free from the 
solar modulation, while it is still sensitive to the local 
magnetic field structure with a spatial scale comparable 
to or larger than the Larmor radius of GCR particles in 
the Local Interstellar Magnetic Field (LISMF). The sky- 
map of the directional anisotropy reported by the Tibet 
AS experiment clearly shows the global feature observed 
with a statistical significance of more than ten times the 
statistical error (TJ. Fig. 1(a) shows the observed GCR 
intensity in 5° x 5° pixels in a color-coded format as a 
function of the right ascension (a) on the horizontal axis 
and the declination (<5) on the vertical axis. The global 
feature of the anisotropy in Fig. 1(a) is successfully 
modeled by a combination of the uni-directional flow 
(UDF) and bi-directional flow (BDF) as 



i: 



GA 
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+ oi|| cosx2(", m : a 2 ,8 2 ) 
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(1) 



where I^m is the GCR intensity in n-th right ascen- 
sion and m-th declination pixel representing the global 
anisotropy (GA), ai± and a^i are amplitudes of UDFs 
perpendicular and parallel to the BDF, respectively. a 2 \\ 
is the amplitude of the BDF, (ax, Si) and (0:2,^2) are 
respectively right ascensions and declinations of the 
reference axes of the perpendicular UDF and BDF and 
Xi (X2) is the angle of the center of (n, m) pixel 
measured from the reference axis of the perpendicular 
UDF (BDF) 0"). This model anisotropy best-fitting to 
the data in Fig. 1(a) is displayed in Fig. 1(b), while 
the residual anisotropy remaining after the subtraction 
of I^m from the data is shown in Fig. 1(c). It is 
clear that this model successfully reproduces the global 
feature in the data. The residual anisotropy in Fig. 1(c) 
is more local and cannot be reproduced by this simple 
model for the global feature fl]. The best-fit parameters 
(ai±, 02||) a U $ti a 2, 82) are listed in Table I. First, 
we note that the reference axis (a 2 ,5 2 ) of the BDF is 
almost parallel to the Galactic plane and very close to the 
LISMF orientation (a = 300.9° and 6 = 32.2°) reported 
by Frisch et al. 0. This suggests the BDF streaming 
along the magnetic field. Second, we note that au_ is 
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Fig. 1. 2D-sky maps of the observed and reproduced GCR intensity. 
Each panel displays the normalized GCR intensity or significance in 
5° X 5° pixels in a color-coded format as a function of the right 
ascension on the horizontal axis and the declination on the vertical 
axis. In this figure, the average intensity in each declination belt is 
normalized to unity. These sky-maps cover 360° of the right ascension 
but cover only 90° of the declination due to the event selection 
criterion limiting zenith angles to < 45° [TJ. The data in 16 pixels 
containing the known and possible gamma ray sources are excluded 
from the best-fit calculation and indicated by black pixels. The white 
curve indicates the galactic plane, while the black curve displays the 
HDP plane suggested by Gurnett et al. 0, which is calculated as 
a plane normal to the orientation of a = 332.1° and S = 35.5°. 
In each panel, the heliotail direction (a = 75.9° and 6 = 17.4°) 
is indicated by a white solid circle on the HDP plane. Each panel 
displays, (a): the observed intensity (-TfJ 6 ^), (b): a component best- 
fit model anisotropy (Z?™) reproducing the global anisotropy (GA), 
(c): the significance of the residual anisotropy remaining after the 
subtraction of 1^1^ from I^J^. In the panel (b), the open triangle 
with an attached character 'F' indicates the LISMF orientation by 
Frisch (a = 300.9° and S = 32.2°) 0, while the open diamond 
with 'B' indicates the orientation of the best-fit BDF (see Table I). 

TABLE I 

Best-fit parameters in I®-^. Note that <5 2 in the upper 

TABLE IS DERIVED ACCORDING TO OUR DEFINITION THAT THE 
REFERENCE AXES AND (02,62) ARE PERPENDICULAR TO 

EACH OTHER |[6). 



H±{%) cn[[(%) a 2 ||(%) 




«(") 


M°) 


0.141 0.006 0.140 


37.5 37.5 


102.5 


-28.9 



much larger than a^i indicating the composite UDF 
almost perpendicular to the BDF. If the orientation of the 
deduced BDF represents the orientation of the LISMF, 
the UDF perpendicular to the BDF orientation implies 
a contribution from the streaming perpendicular to the 
LISMF. Adopting the parallel diffusion coefficient of 
re it ~ 10 29 cm 2 /s by together with the Larmor radius 
of rh ~ 0.002 pc for 5 TeV protons in 3 fiG magnetic 
field, we obtain a Bohm factor 77 = Xvi/tl ~ 1500. 
This large factor implies that the perpendicular diffusion 
coefficient in the LISMF is much smaller than k\\ and 
the observed UDF perpendicular to the LISMF is mainly 
due to the drift flow. The deduced amplitude of the 
UDF (aij_) m Table I is as large as 0.14 %. This large 
amplitude, if it is arising from the drift flow, requires 



a small scale structure for producing a large density 
gradient in the local space. The amplitude of the drift 
flow, expressed as a vector product (B x Vn) between 
the LISMF (B) and the spatial gradient (Vn) of GCR 
density (n), is estimated as 



ai_ 



tl — |Vn| 
n 



1 n 
n L 



L ' 



(2) 



Using the observed ai± — 0.141 % and ~ 0.002 pc, 
we get L ~ 1.4 pc. This implies that we need to 
consider the GCR propagation in a very local structure 
surrounding the heliosphere. As an example for such 
local structure, we consider the GCR propagation in the 
Local Interstellar Cloud (LIC). The LIC is an egg-shaped 
cloud with a volume of 93 pc 3 filled with relatively 
warm interstellar gas [4J. The heliosphere locates inside 
the LIC close to its boundary. There is another cloud 
called G cloud approaching to the LIC and the space 
between two clouds seems to be compressed 0. We 
suggested that the UDF and BDF may be expected 
if the GCR density (n) is lower inside the LIC than 
outside, for instance due to the adiabatic expansion of 
the LIC surrounded by the LISMF. The density gradient 
is maintained by an equilibrium between the inward 
diffusion across the LISMF and the adiabatic cooling 
due to the expansion (6). In this case, the BDF is 
expected from the parallel diffusion of GCRs into the 
LIC along the LISMF line connecting the heliosphere 
on its both ends with the region outside the LIC, where 
the GCR density is higher. The UDF perpendicular to 
the LISMF, on the other hand, is expected from the drift 
anisotropy as described above, if the sufficient density 
gradient is maintained by the adiabatic expansion of the 
LIC. In the next section, we present a simple model 
for the transport of GCRs into an expanding LIC and 
discuss the physical parameters required for reproducing 
the perpendicular UDF observed by the Tibet AS array. 

II. GCR TRANSPORT INTO THE LIC AND 
AN ESTIMATION OF K±. 

We study the transport of GCRs into a spherical 
LIC. The spherically symmetric distribution of the GCR 
density in this LIC is governed by the following radial 
transport equation. 



dF 
ds 



HQ 



dx 2 x dx 



where F(x, s) is the GCR density at a dimensionless 
distance x from the center of the LIC and at a dimen- 
sionless time s, while 7 = 2.7 is the spectrum index 
appropriate for high-energy GCRs. We convert the actual 
radial distance r and time t to x and s, respectively, as 



and s = log e — 



R(t) 



(4) 



where R(t) is the radius of the LIC at time t and t c 
is an arbitrary reference time. We assume a self-similar 
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TABLE II 

Diffusion Coefficient and a ratio of coefficient in the interstellar and the heliosphere 





K ± 


(cm 2 /s) 


re II (cm 2 /s) 


Kr(cm 2 /s) 


K ±/ K T K±/ K \\ 


The interstellar (5TeV) 


4 


X 10 24 


10 29 


6 x 10 2b 


0.07 4 X 10" 5 


The heliosphere (50GeV) 


3 


x 10 21 


3 x 10 23 


7 x 10 21 


0.4 10~ 2 



expansion of the LIC with radius R(t) and expansion 
velocity V(r,t) given as 



V(r,t) = 



t 



R(t) = ^ 



(5) 



(6) 



with R c denoting R(t) at t = t c which we set at 
the present time. The first and second terms on the 
right hand side of eq. (3) describe the inward diffusion, 
while the third term represents adiabatic cooling due to 
the expansion. Note that the convection term due to V 
does not appear in this equation, kq is a dimensionless 
parameter depending on the rate of cross-field diffusion 
and defined by the cross-field diffusion coefficient n± 
as 

(7) 



V c R{t) 



where V c is the expansion velocity of the LIC envelope 
at t c . We set R(t) to be 2.8 pc according to the reported 
volume (93 pc 3 ) of the LIC |4|. The solution of eq. 
(3) rapidly reaches an equilibrium due to the balance 
between inward diffusion (causing an increase of F) and 
adiabatic cooling (causing a decrease of F). The steady- 
state solution F stsa( ty is obtained from eq. (3) with the 
left hand side set equal to zero and given as 



steady.) = Bmh(ax) / 
x sinh(a) V 



2 + 7 



(8) 



xsinh(a)' V K o 

We use this solution to estimate ko and k±. We first 
calculate the fractional density gradient lOO(dF/dx)/F 
obtained from F in eq. (8) at x = 1 where 
the heliosphere is currently located. Fig. 2 displays 
!00(dF/dx)/ F as a function of kq. The fractional 
density gradient is related to the amplitude a\±_ of the 
perpendicular UDF as 

1 dF\ _ R{t) 



100 x 



.F dx J x=i 



(9) 



Substituting the observed a\± = 0.141 % together 
with r L = 0.002 pc and R(t) = 2.8 pc, we obtain 
a fractional density gradient of 197.4 %, This gradient 
assumes kq = 0.47 (see eq. (8)) which, in turn, gives 
hi_\_ = K V c R(t) = 4 x 10 24 cm 2 /s, if we assume 
V c ~ 10 km/s comparable to the speed of the helio- 
sphere relative to the local interstellar medium. 

III. Summary and discussion 

The sidereal anisotropy of the multi-TeV GCR inten- 
sity has been observed with an amplitude of ~ 0.1 %. 
If we interprete this amplitude in terms of the UDF due 



X 
o 
o 




Fig. 2. The figure which shows a density gradient of analytic solution 
on the vertical axis as a function of kq on the horizontal axis. Red 
lines indicate a density gradient equal to 197.4 % and kq = 0.47. 



to the parallel diffusion of GCRs in Galaxy, we arrive 
at 



ai\\ 



(10) 



With oiy ~ 0.1 % and L' ~ 100 pc comparable to the 
scale size of the galactic arm, we obtain km = Ayc/3 ~ 
10 29 cm 2 /s which is comparable to the value reported in 
[ 3 1 . This simple interpretation, however, does not hold if 
the observed UDF is perpendicular to the magnetic field. 
According to the best-fit analyses of the 2D sky-map 
of GCR intensity observed by Tibet AS experiment, the 
UDF with an amplitude of 0.14 % is almost perpendicu- 
lar to the BDF along the LISMF in the galactic plane 0. 
In this case, we need to consider the GCR propagation in 
much smaller and more local structures as discussed in I. 
As an example of such structures, we analyzed the GCR 
transport into the expanding LIC in this paper, but we 
do not exclude any other possible local structures with 
scale sizes comparable to or smaller than the LIC. If 
the GCR transport in such local structure actually needs 
to be considered, it also means that we cannot use Kn 
of for interpreting the observed sidereal anisotropy. 
Based on our expanding LIC model, we estimate k± as 
listed in Table II and compare with the value for the 
heliospheric modulation of GCRs with lower energies 
0. In this table, we calculate kt as 

T_LC_ 

3 ' 



kt 



(ID 



Hence we derive a ratio of k±/kt = 0.07 for multi 
TeV GCRs in interstellar space, while this ratio is ~ 0.4 
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for 50 GV GCRs in the heliosphere. These values are 
compatible within a factor of six. On the other hand, 
the use of k\\ of would result in a 1000 times 
smaller ratio of Kx/k\\ in the interstellar space than in 
the heliosphere. This extreme difference tends to suggest 
that the value of Kn ~ 10 29 cm 2 /s [0 makes it difficult 
to interpret the observed sidereal variation in terms of 
GCR transport in very local structures. 
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